gusucode.com > MATLAB与控制系统仿真实践 源代码 > MATLAB与控制系统仿真实践 源代码/第十一章/eg11_2.m

    num0=[1 3];
den0=[2 4 5 8 10];
G=tf(num,den);
Gc=feedback(G,1);
[num,den]=tfdata(Gc,'v');
r=roots(den);
disp('系统闭环极点：');
disp(r)
a=find(real(r)>0);
b=length(a);
if b>0
disp('系统不稳定.');
else disp('系统稳定.');
end